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Abstract. Starting from a classical- mechanics stochastic model encoded in a 
Langevin equation, we derive the natural diffusion equation associated with three 
classes of multi-scale spacetimes (with weighted, ordinary, and 'g-' Poincare 
symmetries). As a consistency check, the same result is obtained by inspecting the 
propagation of a quantum-mechanical particle in a disordered environment. The 
solution of the diffusion equation displays a time-dependent diffusion coefficient and 
represents a probabilistic process, classified according to the statistics of the noise in 
*Z~ ' the Langevin equation. We thus illustrate, also with pictorial aids, how spacetime 

f - ^ ■ geometries can be more completely catalogued not only through their Hausdorff and 

£ — \ spectral dimension, but also by a stochastic process. The spectral dimension of multi- 

fractional spacetimes is then computed and compared with what found in previous 
studies, where a diffusion equation with some open issues was assumed rather than 
derived. These issues are here discussed and solved, pointing in favour of the model 
with g-Poincarc symmetries. 
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1. Introduction 

The concept of spectral (or fracton) dimension d$ of spacetime descends from the 
analogous one for a set, be it a fractal or a Riemannian manifold [U El EJ HI El [6]. 
The idea is to have a test particle diffuse in the set (a non-relativistic 'space'), and see 
how it behaves. This process is governed by a diffusion equation, which for a Brownian 
motion reads 

(d t -K 1 V 2 x )P(x,x',t) = 0, P{x,x',0) = 6{x-x'), (1) 

where the first-order derivative d t is the diffusion operator in time t, Ki is a constant, 
V^ (called spatial generator in probability theory) is the natural Laplacian in the given 
space (eventually containing metric structure) acting on the x dependence of the solution 
P, and x' is the initial point where the probe, a pointwise particle represented by the 5 
initial condition at t — 0, starts diffusing. The process has a probabilistic interpretation 
when P is a probability density function (PDF), P > for all x. In particular, the 
random position field X(t) associated with equation (pQ) is a Brownian motion. The 
scaling of the variance (i.e., how the mean squared displacement (X 2 ) of the process X 
increases in time) is related to dg. 

To make sense of this picture in a spacetime context with D topological dimensions 
(one timelike and D — 1 spacelike), one must make a few strong assumptions on this 
diffusive process: 
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(A) Replace the Laplacian V 2 with the covariant Laplace-Beltrami (or d'Alembertian) 
operator □ = V^V^ 1 in imaginary time. Thus, there are D directions parametrized 
by coordinates x^, \i = 1, . . . , D, where x D = it and t = x° is time in Lorentzian 
signature. In order to keep the probabilistic interpretation of the diffusion equation, 
the spatial generator is assumed to be an elliptic operator, hence the requirement 
of Euclideanization. 

(B) Time is thus treated on equal footing with the other coordinates, in conformity 
with the spirit of general relativity. Consequently, the variable t in the diffusion 
equation (CQ) is replaced by an abstract evolution parameter a everywhere. 

(C) In particular, a has the same dimensionality as t (it is a length or time scale, 
[a] = — 1 in momentum units), and the diffusion operator d t in ([1]) is simply replaced 
by d a . 

This construction may be unsatisfactory for various reasons. In general-relativistic 
systems, the diffusion equation is expected to have a covariant form, which should 
survive even under the simple replacement procedure (B)-(C). Yet, a non-relativistic 
diffusion equation is typically assumed for curved manifolds. Second, statistical- 
mechanics time t is Euclideanized in the spatial generator □, but not in the diffusion 
operator (which would result in a Schrodinger equation). The general feeling is that, 
even when the diffusion equation is more or less well motivated in the context of 
statistical mechanics, it carries a considerable level of arbitrariness when promoted to 
the diffusion equation associated with a given spacetime. In this sense, the diffusion 
equation (and, hence, the spectral dimension d$) of a spacetime cannot be truly derived 
from solid first principles. 

For ordinary spacetimes, the proof of the pudding is in the eating, and the results 
of the procedure (A)-(C) are reasonable enough not to require a revision. We can 
better appreciate the problem, however, when moving to multi-scale spacetimes, in 
particular those with multi-fractional measure [71[8ll^[in|[TT|[T2|[T3|[TI|[T5|[T6]. 
These spacetimes [7J [17] have been recently introduced as realizations of geometries 
with anomalous (or even fractal-like) properties [T8| [T9| [20] . which appear in certain 
regimes of many quantum-gravity models (a canonical list of references can be found, 
e.g., in [91 [131 CHI El [22] ). The spectral dimension of multi-fractional spacetimes has 
been computed in [191 El HE] using a certain Ansatz for the diffusion equation, and found 
to be anomalous (ds ^ D). Such diffusion equation, however, was simply assumed and it 
turns out to have some issues, which we discuss here. We completely revisit the problem 
(and address these issues) from the more fundamental point of view of statistical and 
quantum mechanics. Although the resulting diffusion equation will be different from 
the old assumption, the spectral dimension will still be anomalous, i.e., different from 
the integer topological dimension of spacetime. However, one must impose certain 
conditions in order for these anomalous spacetimes to be fractal in the usual sense. In 
practice, only one multi-scale model (with so-called g-Poincare symmetries) will satisfy 
all these conditions. 
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As in [121 EH [23], one of the key messages of the present work is that the spectral 
dimension is only one of the many possible ways to characterize a spacetime geometry: 
it is well known that an anomalous correlation function, whose exponent is governed by 
ds, can be obtained by a number of inequivalent diffusion processes [21]. To remove this 
degeneracy, it is important to gain insight into the type of stochastic process underlying 
the diffusion equation. Thus, spacetime geometries can be classified not only according to 
their Hausdorff and spectral dimension, but also by a stochastic process. More generally, 
there exists a well-furbished 'alternative toolbox' [3 [131 03] mutuated from various 
disciplines of mathematics and mathematical physics (chaos theory, probability theory, 
transport and percolation theory, statistical mechanics, complex systems, multi-fractal 
geometry) which is still to be fully tapped into by the quantum gravity community, 
where discussions are often limited only to a few geometric indicators (Hausdorff and 
spectral dimension). We will provide a graphic representation of the possibilities of this 
arsenal by plotting the trajectories of random walkers associated with different fractional 
geometries. 

The plan of the paper is as follows. In section [2] we review the case of ordinary 
Minkowski spacetime, starting from the derivation of the diffusion equation under 
the assumption of having a random walker with a certain stochastic interpretation 
(Brownian motion). This standard setting will introduce the main tools of statistical 
mechanics and probability and transport theory needed also in multi-scale spacetimes, 
the latter reviewed in section [3] In particular, in section 13.21 we analyze the diffusion 
equation previously assumed for these spacetimes and its related problems. Section H] 
presents the new results for various classes of fractional spacetimes, using and extending 
the 'alternative toolbox' introduced in section [2J The spectral dimension of multi-scale 
geometries is computed in section [51 where it is shown that the theory with weighted 
Laplacian can have an ultraviolet regime ds ~ provided one interprets the scale in 
the measure as one signalling a transition between a 'fuzzy' and a continuous regime. 
Section [6] is devoted to discussion. Some of the results related to anomalous scaling can 
be obtained independently (and in agreement with the methods of the main text) from 
quantum mechanics; these digressions are contained in two appendices. 

2. Brownian motion, diffusion equation and spectral dimension 

The spectral dimension of a smooth classical manifold is an indicator of the geometry and 
topology of spacetime. It is obtained by letting a test particle diffuse on the Euclidean 
version of the manifold, and calculating the probability to find the probe again at the 
starting point after some diffusion time a (not to be confused with coordinate time t, 
which is treated on the same ground as a spatial coordinate). This return probability is 
a function Z(o~) of diffusion time; if Z ~ a _ds//2 is a power law, the spectral dimension 
is the exponent ds- 

The test particle is assumed to follow a random walk of Brownian type. This is 
because the associated probability density function yields, in the absence of curvature, 
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the correct spectral dimension ds = D, coinciding with the topological dimension of 
spacetime. The diffusion equation of Brownian motion [25l [26| [27] (see [281 129] for 
extensive presentations) can be derived from a Langevin equation. We review Brownian 
motion and its diffusion equation in a non-relativistic space in sections 12.11 and 12.21 
Promotion to the spacetime picture according to (A)-(C) will take place in section 1231 

2.1. Brownian motion 

Let X(t) be the random variable denoting the position of a particle at time £; in D 
dimensions, it is a vector. A Wiener process (or Brownian motion) X BM (t) with initial 
condition X BM (0) = x' is such that: 

(i) X BM is continuous in t almost surely (i.e., with probability 1); 

(ii) the increments of X BM are uncorrelated, meaning that X BM (t 2 ) — X BM (ti) is 
independent of X BM (t 4 ) — X BU (t 3 ) if the intervals (t 1 ,t 2 ) and (£3,^4) do not overlap, 

(*i,t 2 )n(t 3 ,t4) = 0; 

(iii) X BM is governed by a Gaussian distribution. 

To make the last requirement explicit, let us introduce a dimensionless white noise ri(t), 
that is, a Gaussian random field such that 

(//(*)>, = 0, {7l(t)vW)r, = Kl&(t-t)> ( 2 ) 

where K\ is a constant with engineering dimension [k-\\ = — 1. The noise is called 'white' 
because its spectrum (the Fourier transform of its two-point correlation function) is a 
constant independent of the frequency. Here, angular brackets with subscript r\ denote 
the average over the stochastic background, i.e., with respect to the PDF u(r),t) of the 
process: for any /, (f(r],t)) v := f d D rj u(rj , t) f (rj , t) . In the absence of external forces, 
the Langevin equation of Brownian motion is 

d t X BM (t) = V {t) . (3) 

By definition, it represents a Wiener process as the integral of a white noise: 

X BM (t):=x'+ f dt'r)(t'). (4) 

Jo 

This way, the differential dX BM is well defined even if the trajectory X BU (t) is nowhere 

different iable in the ordinary sense. Two key features of Brownian motion are that it 

obeys the scaling property 

X BM (\t) = \^X BM (t) , (5) 

and that it possesses stationary increments, i.e., 

<[x(t)-x(t')] 2 > = (x 2 (; -;')>• (6) 

Here, angular brackets without subscript denote the average with respect to the PDF 
P of the process, (f(X, t)) := / d D X P(X, x', t)f(X, t). 
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2.2. Diffusion equation 

A Fokker-Planck equation (namely, the diffusion equation) for the PDF of Brownian 
motion can be derived starting from the Langevin equation [281 129] 

mX + m 7 X + U'{X) =F, (7) 

where [m] = 1 = [7], [F] = 2, and dots and primes denote derivatives with respect to, 
respectively, time t and X. This equation describes the motion of a particle X(t) of 
mass m in the presence of a friction force (proportional to the mass of the particle, a 
constant coefficient 7 and the velocity V = X), an external potential U, and a random 
force F, representing the pushing around of the particle by the medium. A writing of 
equation ([7]) where the left- and right-hand sides are dimensionless is 

Lx + x + u ^l = v , ,:=-L. (8) 

7 7717 7717 

There are various methods to find the diffusion equation. One takes the Langevin 
equation (jHJ) exactly and solves it. Consider for instance the case without potential, 
U = 0. Integrating twice the Langevin equation V + '-/V = 777 starting from time t — 0, 
one obtains 

X(t) = f dt'V(t') 
Jo 

= x' + V oX (t) + 7 f dt" x(t - *>(*") , (9) 



where Vb = V(0) is the initial velocity and 

I _ e -7(*-t") 

X (t-t"):= . (10) 

7 

Expression (Q holds true also in more general situations with a more complicated 
function x, for instance when U ^ 0, when the friction term rwyX in ((7|) is replaced 
by a non-local operator and the process is non-Markovian, and for general colored 
noise 77. If the stochastic properties of rj are known, the PDF's resulting from these 
generalized Langevin equations (GLE) [23 EDI ED E21 ESI El E3 ESI EZ! can be found 
by elegant techniques, for instance using the characteristic (or generating) functional 
approach [28], [HH [33] , or the functional derivative approach [32], or the conservation- law 
approach [29]. In particular, the diffusion equation for Brownian motion (also known 
as the Smoluchowski equation when U = 0) is obtained from equation (Q under three 
assumptions: (i) that 77 is a Gaussian white noise, equation ([2]) (thus, also X and the 
velocity field V are Gaussian); (ii) in the limit Vq — > (Maxwell distribution of velocites); 
(iii) assuming to look at times much larger than the relaxation time t rt = I/7. The latter 
corresponds to taking \ ~ const from the start. Namely, setting U = and dropping 
the second-derivative term in (jSJ), we have X w 77, which is nothing but ()3]). 

In fact, for our purposes it suffices to apply the conservation-law approach [29], valid 
for Markovian systems, and to make direct use of equation ([3]). Although the following 
steps are well known, we report them to illustrate the relation between Langevin and 
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diffusion equation for later use. The starting point is the conservation of the probability 

density, 

(1) = f d D XP(X,x',t) = 1 Vt. (11) 

Differentiating this expression with respect to t implies that dtP is proportional to a 
total divergence, in particular (as it happens in fluid mechanics, statistical mechanics 
or electrodynamics) to the divergence of its flux X P. Thus, the conservation law is 
d t P + cVx ' (X P) = 0, for some constant c. From equation (J3J), we can rewrite this 
expression as 

d t P + cV x -(vP) = 0. (12) 

Integrating in time, one has 

P(X, x', t) = P Q (X, x')-c f dt' V x ■ W) P(X, x', t')] , (13) 

Jo 

where Pq(X, x') = P(X, x' , 0) is the initial condition. For times prior to t, the functional 

P depends implicitly on the noise r\. Plugging equation (fT3|) back into f lT2|) . we find 

0. (14) 



d t P = -cV x -(vPo) + c 2 V 



V f dt'Vx ■ {r]P) 

Jo 



We take the average over the noise and call P the stochastic average of the PDF (which 
depends on r], via equation (Ti~2l): 

P(x,x',t):=(P) v . (15) 

The Gaussian statistics of the white noise now leads to the desired result. In fact, 
equation (J2J) implies that the average of the first term in the right-hand side of (1141) 
vanishes (X and rj are independent variables), while the bracket in the second term of 
( TT4l) is made of two contributions: 

f dt'Vx ■ [{v(tW)) v P(x,x',t')} + I* df[ V {f){TiW x ■ P(X,x\f)) v ] , (16) 

Jo Jo 

since P does include noise factors at time t'. The first correlation function gives a 
S(t — t'), the second a 5{t — t") for t" < t'. But t > t" strictly, so the second contribution 
vanishes and, setting c = 1 without loss of generality, we are left with 

(9 t -/t 1 V')P(x,x / ,t) = 0, (17) 

where we have thrown away higher-order noise terms so that the Laplacian acts on the 
spatial coordinates x. 

The diffusion coefficient K\ is measured in units of (length) 2 /(time). Often it is 
effectively absorbed in the definition of a length variable 

i(t):= Visit. (18) 

In general, the diffusion coefficient is determined by the stochastic process underlying the 
diffusion equation (e.g., [281 EH1 [39]). Having averaged over noise, the average (•) with 
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respect to P and the one with respect to P coincide, (f(x, t)) := f d D xP(x, x', t)f(x, t), 
and they shall be denoted in the same way without subscript. In particular, one can 
show that the mean squared displacement (or second moment, or variance) from x' = 
is, in D dimensions, 

(X 2 BM (t)) = 2Dn 1 t. (19) 

The Laplacian V 2 in equation (TT71) acts on the x dependence of the solution P, 
but in this particular case the latter depends only on the distance r = \x — x'\ between 
the initial and final point. In ordinary Euclidean space, the solution is the Gaussian 
distribution 



P(x, x', t) = ui(x, x 1 , t) := — -g , (20) 

[4:7ll 2 {t)\ 2 

where \x — x'\ 2 = (xi — x[) 2 + . . . + (xd — x' D ) 2 and the squared length £ 2 is the 
dispersion of the Gaussian. Consistently with ( ITTT) . the solution is a probability density 
with normalization 

d D xu 1 (x,x',t) = 1, (21) 

■<x> 

meaning that at all times or scales the pointwise test particle (represented by the delta 
initial condition ui(x,x',0) = 5(x — x')) can be almost surely found somewhere. 

2.3. Spectral and walk dimension 

So far we have interpreted the variable t, which has dimension [t] = —1, as a time 

parameter. Moving away from transport theory, in the context of spacetime theories 

t — > a is simply a length, representing the characteristic length scale I (equation ([TBI 

with t — > a) at which one is probing the geometry. The Laplacian is then replaced by 

the curved Laplace-Beltrami operator in Wick-rotated spacetime. 

When a stochastic process is staged on a smooth manifold with boundary, the 

return probability (or heat kernel, or partition function) Z(a) is simply defined as the 

functional trace of the solution P, i.e., its integral over the volume of the set when 

x = x'. Then, the first term in the Seeley-DeWitt expansion (e.g., jlQj HU H2]) is of 

the form Z(a) = V /(A-KKiO~)~ ds l 2 + . . ., where V is the volume of the set. In the present 

case, however, this volume is infinite and it is customary to define the return probability 

as the trace per unit volume: 

, N Z(o) fd D xP(x,x,a) , . 

PW V- ^ D y> . (22) 

This definition coincides with the spatial average of the return probability density 
P(x,x,o~). Since P(x,x,o~) = Ux(x,x,o~) is constant in x for a Brownian motion, the 
denominator in equation (122]) exactly cancels the divergence in the numerator and one 
ends up with V(cr) = [Air£ 2 (a)}~ D ^ 2 oc a~ D ' 2 . The spectral dimension is simply d$ = D 
and it determines the decaying law of the return probability. In general, it is defined as 

n d\nV(a) 

ds(a) = - 2 ^hT^' (23) 
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or, more formally (due to a hidden divergence), 

, . . d In Z(a) , . 

w— 2 -^- _ (24) 

Equation (124"1) is equivalent to the alternative definition ds : = — dlnZ(£)/dln£ only for 
normal diffusion, but in general the correct one is (124)) |1] 121 IH [6] . Dimensionally, 
the mean squared displacement (X 2 (a)) is always proportional to some squared length 
£ 2 (a). 

These definitions are meaningful in a spacetime geometry context provided o (or 
£) is interpreted as a measured length scale. Then, d$ is the scaling law of the return 
probability when a test particle is let diffuse on the manifold. In this respect, one 
should notice that curvature and topology effects do modify the value ds = D at scales 
a larger than the characteristic curvature radius. The classic example is the sphere: 
at o ~ 0, the probe locally feels a Euclidean plane, but as time passes the compact 
topology 'helps' the particle going back to the initial point, so that ds ~ in the limit 
a — > +oo. Therefore, when the spectral dimension of spacetime is computed, the tacit 
understanding is that one is looking for the local geometric properties of the manifold, 
so that in all the cases where Z is not expected to be a simple power law due to topology 
or curvature, it is more convenient to define 

d s := -2 lim ^ . (25) 

o-->o+ In a 



Clearly, when Z is a power law this expression coincides with fl24[) and ds is constant. 

In more exotic scenarios of quantum gravity the effective diffusion equation can 
be very complicated [13] , but at least in semi-classical regimes it must admit a positive 
semi-definite solution P. Otherwise, the operational definition of the spectral dimension 
(I24p as the scaling of the return probability associated with a random walk would be 
lost. The requirement P > was advocated, in particular, in [T3~| [23]. 

The variance defines the so-called walk dimension: 

(X 2 ) oc a 2/dw . (26) 

From equation ( 1T9|) . it follows that <i\y = 2 (normal diffusion) for ordinary Minkowski 
spacetime. In general, the spectral dimension ds precisely determines the anomalous 
scaling in a, since there is a relation between Hausdorff, spectral, and walk dimension 

d w = 2^. (27) 

Consistently, for Minkowski spacetime du = D = ds- This relation holds for all fractals 
and can be understood from (1261) in various ways [6]. A particularly simple one is to 
consider two random-walk systems characterized, respectively, by a length £ (the root- 
mean-square displacement £ ~ -\/{X 2 ) of the walker) and a rescaled one £ — > X£. The 
energy density of states pe(E) is an extensive quantity proportional to the Hausdorff 
volume of the system, pe(E) ~ £ dli . Therefore, under rescaling £ — > X£ one has 

Pxe (E) = \ d «p £ (E) . (28) 
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On the other hand, from the Schrodinger equation one sees that energy E is 
dimensionally conjugate to time, but from equation (126]) (er ~ £ dw ) there follows that 

E u = \- d ™E e . (29) 

For probability to be conserved under rescaling, the density of states in the original 
and rescaled system must be related to each other by p\£(E\e)dE\ E = p^E^dE^, which 
implies, from equation ([29"]) . 

pxi(E xe ) = \ d ™ Pi {ES) . (30) 

The solution of equations (I28p -( 13"U]) is (removing the subscript £ from now on) 

p(E) ~ £<V<iw-i =: E ds/2 - x , (31) 

where in the last step we defined the spectral dimension in analogy with the ordinary 
Euclidean-space expression p(E) ~ E D I 2 ~ X . Via quantum mechanics, one can then show 
(not so straightforwardly) that the definition in ( 1HT|) of ds coincides with the one from 
the diffusion equation. 

3. Multi-scale space-times 

A generic spacetime M. with scale-dependent geometry [18] [7] can be defined by 
equipping ordinary Minkowski spacetime with a measure dg(x), replacing the ordinary 
Lebesgue measure d D x. A non-trivial metric structure independent of the measure can 
be added to describe manifolds with curvature, but we shall ignore it here, as we are 
interested in how a change of the differential structure (and of momentum space) affects 
the properties of spacetime. For technical reasons [El [16] , it is convenient to concentrate 
on factorizable measures, 

D-l 

dg(x) := d D xv(x) = dtv Q (t) dxw(x) := dtv (t) ] J dx i v^x 1 ) , (32) 

i=i 

where the D functions t> M can be all different. By definition, coordinates have dimension 
of lengths ([x] — — 1 in momentum units). A further assumption is that the measure 
weight be positive semi-definite, v^ = v^x^) > for all p = 0, 1, 
of factorizable measure is the fractional one, 



..,D-1. I 


^.n example 


l^lQ^-l 


(33) 


r K) 



dg a (x) = d D xv a (x), v a (x) = ]Jt' a (x M ) := ]J 

where T is the gamma function and a M are D real parameters ('fractional charges') in 
the range < a M < 1. In the simplest 'isotropic' case, a^ = a for all p. In general, we 
will call a the average fractional charge 

a:=— J]a M . (34) 

Fractional measures are not unique (the position of the singularity parametrizes an 
infinite class, there exist non- factorizable versions as in the Riesz integral, and so on [5]) 
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but they epitomize the simplest class of continuous measures with anomalous scaling 
law Q a {\x) = X diI g a (x), where ofa is the Hausdorff dimension of spacetime. For equation 
( 13"3"j) . for instance, 

d H = ^a fi = Da< J D. (35) 

To obtain a scale-dependent geometry, it is sufficient to sum over a finite set of charges 
a [7J [H] • Keeping factorizability, the multi-scale measure weight is 



u.(x):=]Ju,(a^):=] 



5^0nUa n M 



(36) 



where g n > are dimensionful coefficients depending on a hierarchy of length scales 
£ n > 0. For example, an isotropic spacetime with Hausdorff dimension du ~ D in the 
infrared (IR) and g?h ~ -D«* < D in the ultraviolet (UV) is characterized by two charges 
< «i = a* < 1 and «2 = 1 and one fundamental length £i = £* defining large and 
small scales. The measure weight then reads, for each coordinate, 

v,(x) = l + £l~ a *v a ,(x). (37) 

The volume of a D-ball of radius R can be easily calculated and is of the form 
V (D) ~ t?p Dtl (R/£*) D + n D ^(R/Q Da *}, where fi D , Q is the fractional volume of a 
unit ball (£Id,o> = ^D,i/[r(a* + 1)} D if centered at the origin) and we have thrown away 
off-diagonal terms in the integration. When R 3> £*, the first term dominates and the 
ball volume scales as usual; otherwise, the second term dominates and the UV scaling 
is anomalous. In particular, if a* = 2/D spacetime has du = 2 in the UV. 

A more fundamental version of fractional spacetimes is endowed with complex 
fractional measures [TJ [91 [11]. By a suitable choice of coefficients, one can construct 
measures with log-oscillations, of which the present real-order measures are nothing but 
the zero mode. Log-oscillating measures describe spacetimes with discrete symmetries 
and a hierarchy of ultra-microscopic characteristic scales, smaller than those appearing 
in the real-order multi-fractional measure (136]) . At the bottom of this tower, these 
is a fundamental length possibly identified with Planck's length [11]. The interest 
in these models lies in the phenomenological applications of such a scale hierarchy in 
quantum gravity, as well as on the fact that complex fractional measures better represent 
deterministic fractals jl3]. When defining the spectral dimension for these geometries, 
it is necessary to average over the log-period of the measure in order to get a meaningful 
observable [9]. The only surviving contribution is the zero mode, i.e., a power-law 
with real fractional charge. Therefore, all information about the spectral dimension is 
enconded in the models considered in the present work, and we do not need to include 
log-oscillations in the discussion. 

3.1. Laplacians 

We consider three classes of multi-fractional spacetimes, each characterized by a different 
symmetry for the Lagrangian density and, hence, by a different Laplace-Beltrami 
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(Laplacian in short) operator. In the remainder of this section, we leave the form of the 
weight v(x) arbitrary while assuming the factorizability and positivity conditions. 

3.1.1. Weighted Poincare symmetries With respect to the natural scalar product on 
Ai, one can construct a class of self- adjoint operators KL vn of order 27 generalizing the 
d'Alembertian [13]. In particular, the second-order operator belonging to this class is 

m 

1 



K v := rTV^V,,, V^ := -== d M y/v(x) ■ , (38) 

y/v{x) L J 

where rj^ = diag(— , +,•••,+) is the Minkowski metric. A field action endowed with 
this operator is invariant under deformed Poincare symmetries, where the generator 
algebra is the usual one in the free case, but the algebra elements do not generate the 
standard Poincare transformations [16]. In the interacting case, the algebra itself is 
deformed. 

The eigenfunctions of K. v are (here k 2 := k^ = —k 2 + ^ i= ^ kf) 

1 Jik-x 

e(k,x) = — = 5- > JC v <e(k, x) = — k 2 e(k, x) , (39) 



'w(k)v(x) (2vr)- 

where w is the measure weight of momentum space. The normalization is chosen so 
that one can write an invertible momentum transform: 



-00 



f(k):= d D xv{x)f{x)e*{k,x), (40) 

J— 00 

r+00 

f(x)= / d D kw(k)f(k)e(k,x). (41) 



The multi-scale generalization of the Dirac distribution is then 

5{x-x') f D 



5 v {x, x') := / ' = / d D k w{k) e*{k, x)e{k, x') , (42) 

y/v(x)v(x') J 

with a similar expression for 5 w (k,k'). The operator ( |38l) was considered in [81 191 fTU| 

m m m ca m- 

3.1.2. Ordinary Poincare symmetries In another scenario [TBI [T9l [20l [8], [9] the 
Lagrangian density possesses ordinary Poincare symmetries and the Laplace-Beltrami 
operator is simply 

D:=7r0A> ( 43 ) 

which is not self-adjoint with respect to the scalar product with measure weight v(x). 
In fact, 

□ f = IC v :=^-n[v(x)-] . (44) 

Hence, a (for instance) scalar field model with quadratic-like kinetic term —d^<f)d fl <j)/2 is 
inequivalent to the one with Gaussian-like term 0D0/2. In this case, the field theory has 
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some ordering prescription and there is no direct definition of a self-adjoint momentum 
operator (and there may appear complications with the micro-causality structure of the 
theory [5]). Correspondingly, it is not clear whether an invertible momentum transform 
exists. Still, the operators ( 143]) and (144]) will be of interest in what follows. 

3.1.3. q-Poincare symmetries The action measure can be recast as dp = d D q, where 
dg M = dx^v^{x^). Notice that this writing, as well as the most general one q(x) = g(x), 
is equivalent to equation (133]) only in the sense of distributions. In the case of a fractional 
spacetime with fixed dimensionality, for each direction one has 

q»(x*) = Q a {x») = 6 " • 45 

r(a^ + 1) 

Obviously, the Hausdorff dimension is not affected by a change of variables [8] , and 
it is equation (|35|) . du = Da. This happens because by definition the momentum space 
of the theory is conjugate to position space in x, not in q. Consequently, the units of 
the x coordinates are [x M ] = — 1 while the g's have anomalous scaling, \q^\ = —a^. This 
feature guarantees that the theory is not trivial [15] even if one defines it such that in 
q variables it is formally identical to the ordinary one, including the Laplace-Beltrami 
operator 



d d 1 

□ (x) = rf v = rf v d 

dq^(x) dq v {x) v^(x) 



-d v 



(46) 



y v \x) 

which is, therefore, self-adjoint under suitable boundary conditions. In the last equation 
we omitted the indices in the arguments of q and v . 

The measure, the Lagrangian, and the action as a whole possess 'g-Poincare' 
symmetries, i.e., they are Lorentz and translation invariant under transformations over 
the g's, which are then non-linear transformations over the x's (SUH]: 

g"V) = A M ,<f (x) + a? . (47) 

3.2. Previous Ansatz for the diffusion equation 

The explicit calculation of the spectral dimension d§ was performed in jH] for fixed 
dimensionality and in [9 _[T3"j IT6] for multi-fractional spacetimes, in the case of the 
weighted Laplacian fl35]) i|l We recall here the result for an integer-order diffusion 
equation and an isotropic fractional measure (a^ = a, fixed dimensionality). The multi- 
scale anisotropic case is more complicated and adds nothing to the main point we wish 
to make here. 

The Ansatz adopted for the diffusion equation was based on two assumptions: 
(a) since a is a fictitious variable, the diffusion operator d a should not reflect the 
differential structure of fractional spacetime and can be left unchanged; (b) as the 
spatial generator in the ordinary diffusion equation ( IT71) coincides with the self-adjoint 

| A heuristic estimate of ds for the ordinary Laplacian can be found in JT9] . 



CONTENTS 14 

Laplacian appearing in the action, so should it be in the fractional case. This singles 
out the diffusion equation 

{d -K 1 K v )P{x,^,&) = 0, P(x,x',0) = 8 v (x,x'). (48) 



The first problem arises because the solution P(x,x',a) = U\{x, x' , a)j yv (x)v(x') is 



not normalized to 1. It is easy to see that P := ^Jv(x')/v(x)P = U\/v(x) obeys the 
same initial condition, is normalized to 1, and is solution of 

(d a -K 1 )C v )P(x,x',a) = 0, (49) 

rather than of equation fj4*8|) [T3] . The operator K v is not self-adjoint, which requires 
to abandon assumption (b). This is not an issue, since Fokker-Planck equations are in 
general not self-adjoint [29]. One should now explain, however, why the spatial generator 
carries weights v instead of y/v, but for the time being we just ignore this point and 
move on. 

In analogy with equation ( 1221) . we maintain the definition of return probability as 
the functional trace of the solution P per unit volume: 

V(£) = ^- = ^-J dg(x) P(x, x, I) , V H := J dg(x) = J d D x v(x) . (50) 

In fractional spacetimes, however, the spacetime integral in the denominator does not 
cancel the one in the numerator. In fact, the return probability reads 



Id 



D 



X 



V{£) = -^ — w . (51) 

$d D xv(x)^ n p)% 

The ratio of the two spatial integrals is a divergent, dimensionful constant. To get rid 
of it, one must first extract all the dimensional dependence. This is done by defining 
the dimensionless coordinates x M = x^/i, so that V(£) = C£~ Da . The constant C is 
formally divergent but it can be regularized [13] and, since it is dimensionless, be thrown 
away. It is completely immaterial in the definition (|23|) . which would give 

d s = Da. (52) 

Unfortunately, there is an issue with this procedure, which we can state in several 
ways. One is to notice that the regularization trick would introduce a degree of 
arbitrariness in the return probability which, if the underlying stochastic process is 
well defined, should not be there. In other words, the definitions ( 123]) and (|24l) are 
inequivalent in this case, while they should always agree because the overall volume 
prefactor in V does not depend on diffusion time. This would suggest to use (1241) 
or, which is the same, to define V as V = Z/(fd D x) instead of equation (150]) . thus 
modifying the result of the spectral dimension to <ig = D. This is not even a fractal, 
since d$ = D > d^ (for fractals, d$ < du)- The interpretation of V as the spatial average 
of P would also be lost. On the other hand, one might waive away worries that the 
overall constant in the return probability (151]) is regularization-dependent on account 
of two observations. First, in scenarios with fixed dimensionality the regularizion only 
affects the normalization of the PDF, not its signature. Second, when d$ is multi-scale, 
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transient regimes between asymptotic plateaux may also depend on the regularization, 
which is a known feature of multi-scale systems including field theories of quantum 
gravity [T3] . 

Another way to see the problem is to recall that the heat kernel for fractals is of 
the form 

where Vh = L dn is the spectral volume of the set, L is the characteristic spectral length 
determined by the Laplacian (roughly, by the periodicity of its eigenfunctions) and 
dn is the Hausdorff dimension of the set [HI 05]. On the other hand, the divergent 
contribution from the heat kernel Z in (jSTl) is the ordinary integer volume, not the 
Hausdorff volume. This version of the problem can be mildened by recalling that the 
first Seeley-DeWitt coefficient in the heat kernel is quite generically equal to the spectral 
volume but, to the best of our knowledge, it has not been proven as a strict rule for 
non-smooth sets. It is not obvious, however, why this should not be the case also for 
fractional spacetimes, which are a hybrid between continuum geometries and genuine 
fractals. 

A third manifestation of the problem becomes apparent when one identifies its 
origin in the normalization of plane waves ( |39l) : the phases e(k,x) are normalized per 
unit integer volume, not per fractional volume. This is apparent from equation (I42p . In 
yet other words, the number of states is not the fractional volume in equation ( 150]) but, 
rather, 5(0) /v (x). In fact, define the density of states of a system with Hamiltonian 
eigenstates ip k (x) and energy E k as 

w(E) p(E) := I ' d D kw(k)\ij k (x)\ 2 5(E - E k ) , (54) 

where w(E) is the energy measure weight. Then, the total number of states is 

dEw(E)p(E)= [d D kw(k)\Mx)\ 2 = ^r\: (55) 

J v(x) 

where in the last step we used the fact that the states ip k always contain a factor of the 
form 1/ \Jv(x)w(k) (as in the free particle case, where ip k (x) = <e(k,x)). In the bra-ket 
formalism, this is tantamount to starting with the resolutions of the identity 

1* = jAd]^ j ^ D ^{k)\k){k\ , 1* = j^y^Jd D xv(x)\x)(x\ , (56) 

so that |x) = (2n)~ D l 2 J d D x'v (x')\x')(x'\x) and one finds (x'\x) = (2ir) D / 2 5 v (x,x'). 
Since (2 / k) d ' 2 5 v (x, x') = (x'\t k \x), the right-hand side of (|55|) can be recast as (x\x) = 
5(0) /v(x), where 5(0) is the volume of integer momentum space. It is the integer volume 
because the eigenfunctions ip k (x) always cancel the weight factor w(k). Therefore, each 
infinitesimal hypercube d D k (or d D x) contains w(k)d D k/w(k) = d D k localized states 
(respectively, d D x), as in the usual case. 

Ultimately, the problem lies in the diffusion equation itself, equation (1491) . which 
was simply assumed. Starting from the Langevin equation we will get, instead, a 
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different diffusion equation not only derived in closer conformity with the standard 
case of ordinary space, but also capable of overcoming the difficulties outlined above. 

4. Diffusion in fractional space-times 

4-1- Weighted Laplacian 

The natural extension of classical mechanics to multi-scale spacetimes is obtained by 
replacing ordinary time derivatives with weigthed derivatives T> t |14j . For instance, a 
free particle in a quadratic potential has action 



S = dtv (t) 
with equation of motion 



1 IT, \2 1 2 2 

-m[L>fX) moo x 

2 y ' 2 



(57) 



mD 2 x + muj 2 x = . (58) 

Similarly, it is natural to define the Langevin equation for a random variable X in 
fractional or multi-scale space as 

mV 2 t X + m^V t X + U'(X) = F. (59) 

Switching to the spacetime interpretation (A)-(C) spelled in the introduction, one ends 
up with the stochastic equation 

mV\X + mi V a X + U'{X) = F, (60) 

where we switched to diffusion time according to assumption (B). However, now 
assumption (C) is less transparent than in the standard case, since the derivative T> a 
contains a weight v(a) for an abstract evolution parameter and, a priori, it is not obvious 
whether this weight should have the same functional form as the weight Vo(t). For the 
time being, we leave the form of v(o~) unspecified, calling 1 — (3 = [v(a)} its scaling 
dimension. In particular, in the fractional case we set 

v{a) = v p {a) := ^— , > , (61) 

with P not necessarily equal to the fractional charge cto- 
Making now the field redefinition 

Y{a):=yfi{a)X{a), £(a) := — y/^a)F{a) , (62) 

equation ( 160]) reduces to (jSJ) in the presence of a non-autonomous potential W = 
^Jv(a)U and a noise £: 

h^ + B . Y+ ^^n = (. (63) 

7 m7 

The resemblance with equation (jSJ) is only superficial, since measure units differ: the 
effective variable Y has scaling dimension [Y] = —(1 + (3)/2, while [£] = (1 — /3)/2. 
Notice that if X is a random variable, so will Y, since via) is a deterministic function. 
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Also, at the initial time a = one takes the expression Y(0) = lim^o \f v i a ) X{a) 
formally finite, even if v(0) may diverge. 

A natural assumption is that F is a Gaussian white-noise field but in a fractional 
sense, i.e., with correlation {F(a)F(a'))p oc 5 v (a, a'), where the proportionaly coefficient 
has dimension (mass) 4- ^ {(FF) has scaling dimension 4 and [5 V ] — 0). It is more 
instructive, however, to consider the general case 

(F(a)F(a')) F = (m 7 ) 2 /^ a^S^a, a') =► (£(<r)£(<r')>* = k p>v a^5(a - a') (64) 

for all a and a', where v is a real parameter and K v> p is a constant with dimension 
[ K /3,u\ = v—fl — l. Repeating the same calculation as in section [21 one finds an expression 
for the probability density u (we reserve the symbol P for later) satifsying an equation 
of the form (1X41) . with X replaced by Y . Taking the average over the noise and calling 
u(y,y',a) := (u) v , we end up with 

(d a -K /3 , u a u - 1 V 2 v )u(y,y\a) = 0, (65) 



where y := y/v(a)x actually contains a dependence on a. Restoring the coordinates x, 
we finally obtain 

[d (r -K(cr)V x ]ufi tl/ (x,x',(r) = 0, k((j) := Kp <v —— . (66) 

v{a) 

In particular, for the fractional measure (|6T|) k(ct) = Kp^T{[5)a u ~^ . Notice that the 
solution of the diffusion equation (16"6"1) is not the solution of 065|) . since the function 
u[y/v(a)x, y\ a] is neither correctly normalized to 1 nor well-defined for initial points at 
\x'\ < oo. The actual fractional solution up >v {x, x', a) is a Gaussian PDF proportional 
to u(y,y',a): 



g « 2 (ct) 

Uf, tV (x, x', a) = , (67) 

[A7r£ 2 (a)]~ 

where the dispersion £ 2 is 

i 2 (a) := I 2 + ^ da' k(<j') =I 2 + Kp, v f da' ^— , (68) 

and we allowed for an additive constant I 2 . The initial condition up^(x, x 1 , 0) = 5(x — x') 
imposes that £ 2 (0) = 0, which fixed I 2 depending on the measure v(a). In the fractional 
case (l6i~T) . one has that £ = and 

If we had set £ ^ 0, we would have obtained a Gaussian initial spread up tU (x,x',0) = 
exp[— | a; — x'\ 2 /(4:£ 2 )]/(4:ir£ 2 ) D ' 2 of width £, a situation mimicking the idea of a particle 
on a fuzzy manifold with a minimal length [4*6l Wf\ |4*8| |4*9| [50| ISTj . We will come back 
to this point in section [U where a non- vanishing £ is necessary in certain multi-scale 
scenarios. 
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Equation (I66p is not yet in its final form, since we have not discussed the 
normalization of the solution in a multi-scale spacetime. In spacetimes with non- 
trivial measure, the normalization of the PDF is done with respect to the total measure 
weight v(x). The only way to get the result and maintain both the Langevin-equation 
interpretation and the delta initial condition is to define the fractional PDF as 

Pfr„(x,x,a) = , (70) 

v[x) 

which is the solution of 

[d a - K(<j)it v ]Pp >u (x, x\ a) = , Pp,v(x, x', 0) = S v (x, x') . (71) 

Since up >v obeys the self-similarity relation 

u^ {1+V ~ m ^ \ {l+ ^ )/2 x', Act] = \- D{1+l/ -^ /2 u^(x, x', a) , (72) 

P/3 jV is self-similar with scaling law 

iV[A (1+ ^ )/2 x, A^^V, Xa] = \- Da ^-^/ 2 P^(x, x', a) . (73) 

Quantum mechanics motivates the appearance of the operator K, v in the diffusion 
equation (ITT]) independently (appendix [A]). It also determines, independently from the 
Langevin-equation approach, the natural diffusion time ( !68l) (appendix IB1) . 




4- 1.1. Factorizable- spacetime Brownian motions Set first v — 1. From equations ( 1651 ) 
and (162|) . it follows that X is the stochastic process 

(74) 

This is a plain .Brownian Motion, but in a fractional or, in general, Factorizable multi- 
scale Spacetime (FSBM). We label FSBM in this model of multi-scale geometry with 
weighted Laplacian as FSBM-f . Some trajectories for the weight flBTT) are shown in 
figure CQ The raggedness of the curve and the drift from the average decrease when (3 
increases. 

Together with equation ([5]), equations (161]) and (174j) imply the self-similarity 
property 

2-/3 

^fsbm--u(Act) = A 2 X FSBM _ v (a) . (75) 

The mean squared displacement of FSBM is anomalous, 

(X F 2 SBM _„(cr)) = / d D xv(x)Pp^(x,x',o-)x 2 oc£ 2 {a), (76) 

and, for the fractional case £ 2 (cr) oc cr 2_/3 , one has sub-diffusion when < (3 < 1. 
Moreover, increments are uncorrelated but not stationary, since their distribution does 
not depend on the time interval only: 

(LY FSBM _„(cr) - X FSBM _ v (a')} 2 ) ^ {X 2 SBU _ V (<J - a')} . (77) 
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Figure 1. Thick (black) curve: example of trajectory of the process ([71)1 in a fractional 
spacetime with weighted Laplacian, diffusion measure weight (pT|) . f = 1, and /3 = 1/2 
(top panel), /3 = 1 (middle panel, ordinary Brownian motion), and /3 = 3/2 (bottom 
panel). Thin (blue) curve: example of trajectory of ordinary Brownian motion in 
ordinary spacetime, plotted for reference (notice the different scaling of the vertical 
axes). 
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If v 7^ 1, set £ = o"^ -1 )/ 2 ?], where r\ is a Gaussian white noise ([2]). Then, the 
Langevin equation (from fl63|) ) 9 CT F = £ becomes 

5 ff X SBM (a) = <r^ V (a) , (78) 

where we called Y(a) = X SBM (a) the solution. This is nothing but a scaled Brownian 
motion (SBM) [351 ESI [21] , defined as 

X SBM (a) := X BM «) . (79) 

Its self-similarity property is X SBM (A<r) = A^X^^a). Just like FSBM, SBM is 
Markovian, since the scale transformation a — > a v preserves time ordering for v > 



Again, increments are non-stationary, ([X SBU (a) — X SBM (a')] ) = {[X 



o- - 



X BM (a ,u )] 2 ) = (X% M (a u - a' y )) ^ (X s 2 BM (cr - a')). Figure [2] shows the effect of the time 
rescaling: as v decreases, the trajectory becomes smoother. 

Equations fl74j) and f!75|) are replaced by a fractional-spacetime scaled Brownian 
motion (FSSBM in short): 




(80) 



with scaling law 



^fssbm(Act) — A 2 A / " fssbm _^(<t) . (81) 

Its variance is proportional to £ 2 , which for the fractional case (161]) is 

(* F 2 ss B[ >)>^ 1+ ^- (82) 

We notice that, in an ordinary spacetime, equation (166]) with time-dependent 
diffusion coefficient can describe two different stochastic processes . One is SBM. 
Another is fractional Brownian motion (FBM) (52] [53] (see also [301 EH ESI E3 [351 HE E3 
23]). In this case, the GLE features a Weyl fractional integral of order 7 = (2 + v — (3)/2 
and reads d a X FBU (a) = 9 (T [_ 00 / 7 ?7]((j). Scaled and fractional Brownian motion are 
twins [23], i-e., they share exactly the same diffusion equation ([66]) and the mean 
squared displacement. However, they are physically quite different, since FBM neither 
is Markovian nor has stationary increments. From the point of view of fractional spaces, 
there is no twin problem and no ambiguity. In fact, the natural Langevin equation is 



<9 CT [v^(V)X]=£, (83) 

which is neither the one of SBM nor of FBM. Furthermore, we are dealing with 
a Markovian process, as the derivation of the diffusion equation has shown. The 
Markovian stochastic process associated with equations (183]) and (|66l) is unambiguously 
identified with FSSBM, or with FSBM-w if v = 1. 
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Figure 2. Thick (black) curve: example of trajectory of scaled Brownian motion 
(|79p in ordinary spacetime with v = 1/2 (top panel), v = 1 (middle panel, ordinary 
Brownian motion), and f = 3/2 (bottom panel). The same trajectories represent 
the process ([89| in a fractional spacetime with ordinary Laplacian with v = 1 and, 
respectively, j3 = 3/2,1,1/2. Thin (blue) curve: example of trajectory of ordinary 
Brownian motion in ordinary spacetime, plotted for reference (notice the different 
scaling of the vertical axes). The code implements the algorithm presented in equation 
(A2) of [35]. 
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4-1-2. Spectral and walk dimension For a generic measure v(o~), the spectral dimension 



is 

d s (a) 



m _ d\nZ(£)d\n£ 2 (a) 
d In £ 2 d In a 



m D ^f {a) (84) 

din a 

J da K{a') 

In the case (168]) . the heat kernel is Z(cr) oc <t _d ( 1 +^ _ ^V 2 (where the proportionality 
constant is the total integer volume) and 



<k = D(l + u-p) 



(86) 



The spectral dimension depends both on the choice of weight for a and, via the 
parameter v, on the statistics of the underlying stochastic process. We distinguish 
various cases: 

(i) A natural assumption (also supported by the findings in appendix [B]) is that this 
process is the counterpart of Brownian motion, so v = 1 and ds = D{2 — (3). 

(a) If one further assumes that /3 — 1, one recovers the case of the old diffusion 
equation with non-regularized heat kernel Z, section 13.21 where ds = D. 
Diffusion is non-anomalous, {X 2 SSBM (a)) oc a. 

(b) Assuming that a inherits the measure of time coordinate, (3 = «o, one has 
ds — D(2 — a ) > D and super-diffusion. 

(c) Assuming instead that the charge (3 is equal to the average fractional charge 
( 13^1) . (3 = a, ds = D(2 — a) > D and again one has super- diffusion. 

(ii) The evolution parameter a does not have to inherit the measure of time coordinate, 
and there is no compelling reason why (3 should coincide with the fractional charges 
«o or a, as in cases (b) and (c). One could set /3 — 1 from the start, as in the old 
diffusion equation, but reintroduce a dependence 'from the backdoor', via a scaled 
Brownian motion X = Y with v = ex.. Then ds = Da as in the old diffusion 
equation case with regularized V. 

None of these cases corresponds to a fractal in the usual sense, since from equation 
( 182]) one sees that (X 2 SSBM (cr)) oc a ds ^ D ^ a ds ^ dl1 . The origin of the violation of relation 
( )27l) can be understood by revisiting the argument in section 12.31 with the form of the 
density of states found in section [3721 Apart from the replacement p(E) — > w(E)p(E) 
(which, in fact, is only a redefinition of the density and one can just put w — 1), the 
main difference is in the state counting per Hausdorff volume, which is not proportional 
to do(x) but to d D x. Thus, one should replace g?h i n equations ( 128]) and (|3T]) with the 
topological dimension D. Thus, for this model of fractional spacetimes the relation (J27j) 
is replaced by 

d w = 2^-. (87) 

as 
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We conclude with a remark on the various choices (i)-(ii). The diffusion equation 
is derived, or assumed, in a classical-mechanics context. Here, the only information 
about an anomalous dimension is given by the measure in the time direction, while 
the only information about spatial dimensions is the number D — 1 of particles x l (t), 
which is simply the topological dimension of space. From the point of view of the 
diffusion equation, the time direction is an external, arbitrary parameter a. There is an 
intrinsic element of ambiguity in the whole construction which can be removed only by 
a definition of the statistics of the random noise. Depending on which choice is regarded 
as 'fundamental', the final output will be different. The problem in this model (and the 
next) of fractional spaces is that none of the above choices seems to be well motivated 
when looking at the analogous procedure in ordinary Euclidean space. In section H~3l we 
shall consider another model where this conundrum is apparently solved. 

4-2. Ordinary Laplacian 

Equation (17TT) is not self-adjoint, which is a rather common situation for general Fokker- 
Planck equations [22]. However, in a quantum-gravity setting, or in any case where the 
diffusion equation is aimed at the determination of the spectral dimension of spacetime, 
one may wonder whether the physical and geometric consequences of the diffusion 
equation change if one takes its adjointlj In the case of ( jTTl) . /Cj is simply the ordinary 
Laplacian and the adjoint equation is 

[d a - K{a)Vl]P p , v {x, x', a) = , P PtV {x, x', 0) = 8 v (x, x) . (88) 

This diffusion equation actually corresponds to the scenario with standard integer 
Laplacian. In that case, one can show that /3 is the exponent governing the power 
law of the quantum- mechanical time (appendix |B1). 

The associated stochastic process is the scaled Brownian motion 



X(a) = X FRm ,(a) := X nM (a 



l+v-B\ 



(89) 



In fact, in order to construct the GLE associated with this spacetime, one should account 
for both the differential structure (ordinary derivatives) and the fractional generalization 
of the source £. The first information suggests that the left-hand side of the Langevin 
equation should be of the form d a X. The second information yields, in the simplest 
case v — 1, (£ 2 ) oc 5 v (a,a') and £(<r) oc rj(a)/^v(a), where rj is a white noise. The 
result for general v is then d a X = a^ v ~^^ 2 r] for the weight (1611) . which is (I78p with v 
replaced by 1 + v — (3. The argument is somewhat heuristic, since ordinary derivatives 
mix with weights in this model upon integrating by parts, and it is not obvious how a 
complete GLE should look like. Anyway, it yields the diffusion equation (IHHj) and the 

§ The operation of taking the adjoint of the diffusion equation is not of clear interpretation, so even if 
ds differed in mutually adjoint pairs one should not draw the conclusion that there is a problem with 
them. 
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same anomalous scaling ( 1811) of the previous model. A trajectory of the walker is shown 
in figure [2j 

The solution of equation ( IHH1) is again a Gaussian like ( |67l) . 

„'|2" 



P/3,u = C(x', a) exp 



— \x — X 
A£ 2 (a) 



(90) 



but with a non-standard normalization C depending on the initial point x' . This is 
due to the initial condition in (188 p . different with respect to that of equation f[6"6"j) . We 
impose 



l = C(x',a) / d u xv(x)e «^ . 



(91) 



Consider the case of fixed dimensionality, v — v a . The integral f actor izes and can be 
done analytically for each direction ([51], formulae 3.462.1 and 9.240). Inverting the final 
result yields the normalization constant 



C(x',a) = \ 



rK/2) 
r(a M ) 



[2£(a)} a »$ 



a^ 1 

! '2' 



U' 



/Ml 2 



A£ 2 (a) 



(92) 



where $ (also called \F-\_ or M) is Kummer's confluent hypergeometric functions of the 
first kind: 

° (a) n z n 



$(a;fe; z) := ^ 



n=0 



(&)„ n! 



(93) 



where (a) n = T(a + n)/r(a) is the Pochhammer symbol. To check that the initial 
condition in ( 155|) is respected, we recall the asymptotic limit 



Thus, 



$(a;6;2) 



C^', a) 



n 



r(6-a) v 

r(i/2) 
r(a„) 



.r 



//' 



Otu-l 



{47r£ 2 ) D / 2 v a {x') 



(94) 



(95) 



and one recovers the fractional delta distribution S v . 

Next, we calculate the asymptotic limits of the heat kernel 



Z(a) 



d xv a {x) Pp iV {x 1 x, a) 



d xv a (x) C(x, a) . 



(96) 



This expression is not a power law in a, so in general it will not give a constant spectral 
dimension as in the adjoint case. We estimate ds at small diffusion scales, equation 
( 1251) . Using again (l95l) . we get Z(cr) ~ V(4-7r£ 2 ) _D ' 2 , where V = J d D x is the ordinary 
volume. Thus, one recovers equation (156"j) . 

For large a, $ -► 1, C(x, a) ~ n M [ r («/./ 2 )/ r ("M) rM (^)] _1 an d Z(«r) oc 
VH^ _Da (c) oc <j~ Dat > l+u -f 3 )/ 2 ^ and the spectral dimension is a times smaller than the one 
at small scales. This mismatch can be explained as an artifact of adopting a measure 
corresponding to a geometry with fixed dimensionality. Characteristic scales appear 
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either as topological effects (as, for instance, the curvature radius of a sphere or a 
torus) or in dimensional flows determined by an intrinsically multi-scale measure. Here 
there is no foothold where to establish a scale hierarchy, so it is more natural to take 
equation (125)) as the correct definition of ds when characteristic scales are not expected. 
Therefore, the limit a — \ oo of the return probability has no significance here and the 
result for ds can be interpreted as in agreement with the one from the adjoint diffusion 
equation. Further support to this conclusion will be given in section I5.2[ where we shall 
compute ds in the multi-scale scenarios. In the more realistic case where ds changes 
with the scale, the spectral dimension of the two models is in perfect agreement. 

4-3. q-Laplacian 

This scenario has a standard classical mechanics (ordinary derivatives in a certain 
diffusion time r). All the results of section [2] hold with the following changes: 

(i) The random variable X(a) is now regarded as a composite variable Q{r) = Q[X(t)]. 
In the case of fixed dimensionality, for each direction 

sgnPf MllX^MI " 
Q"[X{r)] = Sg [ pPW )l ■ (97) 

I [atp + 1) 

This expression is invertible, since sgn(X) = sgn(Q) and X = sgn(Q)[r(a + 
l)!^!] 1 ^. The parameter r is by itself composite, r = g(cr), in particular 

T = eM = WTT) < 98) 

in the no-scale fractional case, where a has the dimensionality of a length and (3 is 
a constant. 

(ii) The diffusion equation is 

^-E^)7^(^} P = ' M = -0, [« M ] = - 2a„ , (99) 

where, in the general anisotropic case, the diffusion coefficient is different for each 
direction. One can also define the parameter 

W := Mr(flr)] = yWW » &] = -«m- ( 10 °) 

In the isotropic case (a M = a), there is only one diffusion coefficient k and one 



parameter £(a) := y/KT(a) with scaling [£] = —a. Then, the diffusion equation 
(19"9~|) can be written as 

d d 2 

dl 2 !^) ~ ^ [dq(x»)] 2 j P = °' (101) 

This version of the diffusion equation is also valid in the anisotropic case when the 
dimensionful coefficients are absorbed in the definitions of g M , in which case £ 2 = r. 
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(iii) Thus, as a function of the actual diffusion time a (which is the one with respect 
to which one must derive the heat kernel in order to get d$) the random variable 
Q{a) = Q[r(cr)] is a scaled Brownian motion. In the particular case of equation 
(p7|) . in the first orthant one can identify the stochastic process, which we dub 
FSBM-g, as (up to 0(1) constants) 



v 



X FSBM . q (cr) ~ [X BM (^)] 1/a = [X SBM (a)f a . 



(102) 



The raggedness of the ordinary SBM increases with increasing (3 (figure [2]) . On the 
other hand, with respect to an ordinary SBM the raggedness of FSBM-g is greater 
(1/a > 1) and increases for smaller a; the drift from the average is then amplified. 
There is an interplay of the two effects, but the effect from a (power of Q) is greater 
than the one from (3 (scaled time), even when (3 = a. Therefore, this incarnation 
of FSBM is typically more ragged than an ordinary Brownian motion; see figure El 

(iv) Equation ( 12U|) holds with the coordinate replacement x —> q(x), 



e if-i (") 



P(x, x', v) = \\ , (103) 

where the normalization is automatically correct. 



The heat kernel is of the form (1531). which for the fractional case (1981) reads 



Z(o) = Vh ocV H a-^ D/3 , (104) 



so that the first Seeley-DeWitt coefficient does correspond to the Hausdorff volume 
(this fixes a problem unsolved in the other models) and the spectral dimension is 



d s = D/3. 



(105) 



As in the other theory, the spectral dimension only depends on the topological 
dimension of space and on the fractional charge associated with diffusion time. The 
latter stems from a classical mechanical model which knows nothing about the Hausdorff 
dimension of spacetime, since it only sees D motions Q^. Again, a is just a parameter 
from the point of view of spacetime, and (3 does not necessarily have to be equal to the 
time fractional charge «o- To fix this arbitrariness, we can encode the information of the 
Hausdorff dimension in (3 and define it as the average fractional charge, (3 = a = d^/D. 
This way, it results that d$ = du, as in equation (1521 . 

Even for a general (3, the walk dimension is indeed the one for a fractal, since from 
Q ~ r 1 / 2 there follows X ~ a^ 2 ^, and equation ([27]) holds. 
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X(o-) 




100 200 300 400 



500 



100 



X(O-) 50- 




500 



Figure 3. Thick (black) curve: example of trajectory of the process (|102[) for the 

diffusion measure weight (|61[) . with /3 = 1 = a (top panel, ordinary Brownian motion 
in ordinary spacetime) and /3 = 1/2 = a (bottom panel, fractional spacctime with 
q-Laplacian) . Thin (blue) curve: example of trajectory of ordinary Brownian motion 
in ordinary spacetime, plotted for reference (notice the different scaling of the vertical 
axes). 



5. Diffusion in mult i- fractional spacetimes 



Assume that the distribution v(a) takes the same form of the multi-fractional measure 
weight (136]) . 

TV 

v(a) = v,(a) = J29n^ n ~\ (106) 

71=1 

and for the sake of simplicity consider the binomial case gi = ^ _/3 *, < /3i = /3* < 1, 
g 2 = 1, and j3 2 = 1: 

/ \/3*-l 

v*(<r) = 1 + Ij-J , (107) 

where £* is a fundamental length discriminating between infrared (er ^> £*) and 
ultraviolet (<r <^C £*). Also fix v = 1 in the statistics of the random walker and call 

ft/3,1 = Kp- 
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5.1. Weighted Laplacian 

When the spacetime measure weight v(x) encodes a multi-scale geometry, the theory 
with weighted Laplacian does not fare well, although it can be rescued by an interesting 
modification of the physical interpretation. 

The multi-scale version of the process ( 1731) is obtained by using the distribution 
(I106P or (11071) ; we do not plot the resulting trajectories here. From now on we aim to 
compute the spectral dimension. 

Let 1 < /3* < 2. From the asymptotic behaviour of (11071) . one expects that in the 
small-scale limit the first term dominates and d§ ~ D, while in the large-scale limit 
the second term takes the lead and one recovers d$ ~ D(2 — /3*) < D, i.e., a reduction 
of the spectral dimension with the scale. This is confirmed by an explicit calculation. 
Plugging (11071) in equation flBBl) . we obtain 

da' OT I" 1 & / a x /3 *~ 1 " 

1 + (a'/QP*- 1 



<V)=?+/ ; . rTTT^TT -^ + ^F 



'/3*-l'/3*-l' 
where we used formula 3.194.5 of [51] arid F = 2-F1 is the hypergeometric function 



(108) 



F(a , 6;c;z):= £(^J. (109 ) 

n=0 Wn 

For a <C £*, F — » 1 and ^ 2 (cr) — >• £ 2 + Kpa, thus fixing £ = 0. For z — > 00 one must 
analytically continue via formula 9.132.1 of [54"] : 

F(l, 6; 6 + 1; s) = & ! n ^/ (l, 1; 2 - b; j^j + r(6 + l)r(l - b)-^-~ b . (110) 

For 0, > 1, b = l/(& - 1) > 1, so that 

K/30" (a <C £*) 

2 ^M «-<r (aV^* , , forl<&<2. (Ill) 

(a > t 



From equation (1851) . the spectral dimension reads 



*<" = ^- (112) 

From the asymptotic behaviour (II lip , it immediately follows that 

{ D (a < 4) 

d s ~ < , for Kj3* <2, (113) 

\ D(2 -P*)<D (a > £*) 

as announced. 

Let us now consider the case < /3* < 1 (6 < — 1). Again, one expects a decrease 
of the spectral dimension as the scale increases, from ds ~ -D(2 — /?*) > D to ds ~ D. 
Equation (11081) is valid only for /?* > 1, but we can analytically continue it under the 
provision that £ 2 (0) = 0. At large scales, £ 2 {<j) ~ £ 2 + Kpcr, while for a 3> ■£* we use 
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equation flllOp . this time picking both terms. Noting that r(6)T(l — b) = 7r/sin(7r6), 
from equation f 1 1 8 j) we get 



()2l \ <r<£l* 



F + Kpt 



sin 



0.-1 



KfjCF I O 



2-/3, U 



(114) 



implying that 



/?*-! 



K/,4^^-- (115) 

sin 



/9.-1 

This formula is valid also for ft = 1 + 1/&;, k Eli, where I 2 — > —Kp£*. To summarize, 

{Kpa f cr\ , p x 

2-/3, Vt J {<J < * j , for < ft < 1 , (116) 

P + Kpa (a > 4) 

and we get an interchange of the regimes of ( II 13ft : 

ds~< n ' , forO<ft<l, 117 

although the decreasing behaviour remains. 

The overall scenario is radically different from the 'embedding picture' previously 
associated with models with weighted Laplacians. In that case, one started with a D- 
dimensional ambient Minkowski spacetime M D and had a multi-fractional spacetime 
.M* embedded in it. Then, .M* 'filled' M D only at large scales. This is true also in the 
present case, but only for the Hausdorff dimension. On the contrary, there is a reduction 
of the spectral dimension at large scales: 

I. For < ft < 1 (equation (I117p ). it is as if fractional spacetime gushed out of 
the embedding at smaller and smaller scales. Taking, for instance, D = 4 and 
= 1/4, 1/2, 3/4, one would have a dimensional reduction d$ ~ 7, 6, 5 — > d$ ~ 4. 

II. For 1 < ft < 2 (equation (11131) ). it is as if the embedding were depleted of spacetime 
points at larger and larger scales. For example, for D > 5 and f3 = 2 — 4/D, one 
gets ds ~ 4 in the infrared starting from an embedding with 5 or more dimensions. 

Possible modifications may come from taking a more general statistics for the diffusion 
process [y ^ 1, but this is not natural from the point of view of quantum mechanics; see 
the appendices) or from changing the initial condition from a delta to a Gaussian. As 
discussed in section |4"7l] this can be achieved simply by letting the constant £ 2 to be non- 
vanishing by default. The interpretation then would be the one of (IHlllTlllEJEO], name hy 
that the test particle is not pointwise due to the intrinsic 'fuzziness' of a spacetime with 
a minimal length. To get this type of structure in the multi- fractional context, it would 
be quite natural to identify the characteristic scale £* in the binomial measure with the 
scale at which 'fuzziness' effects become important. In turn, this univocally sets 

£ = L. (118) 
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This would 'straighten' dimensional flow. Take, in fact, < /3* < 1. The only but 
crucial modification to the above calculation of ds is that the constant term in equation 
(11141) does not vanish. Call L 2 the content in square brackets. While the large-scale 
limit of the spectral dimension is unaffected, the small-scale one is 



,8.-1 



. (a < Q , s 

d s ~ { . 120 

^ D (a > L) 



leading to 



Such a scenario somewhat merges multi-fractal]]] and effective non- commutative 
geometries, but in a way different from [TT]. In [11], ^-Minkowski spacetime was 
reinterpreted as the ultra-microscopic limit of a log-oscillating fractional measure. On 
the other hand, here we have changed the initial condition of diffusion (but, implicitly, 
also of the Green function) so that to have the scale in the real-order multi-fractional 
measure play the role of a fuzziness length, just as in commutative effective spacetimes 
coming from a non- commutative geometry where operators are evaluated on coherent 
states [lEJ E7J HE]. However, the matching is only at the level of the spectral dimension 
[50] . since our resolution of the identity is not a Gaussian but a fractional delta. We leave 
a further elaboration of this scenario and a fuller assessment of its physical consequences 
for future study. 

5.2. Ordinary Laplacian 

The normalization condition ( 19TT) now features a multi-scale measure weight v(x) = 
v*{x). With the binomial measure (|37|) and in the isotropic case a^ = a*, one has 

T(a,/2) 2 a *£ a *(a) 



2/_M# , nD 



C(x',a) = { [47rf{a)}- + t 



r(«* 



r l-a, 1 (x >lx ) 2 



n* 



1 c\ 1 



2' A£ 2 (a)_ 



(121) 

where £ 2 {<j) has been computed in the previous sub-section. As before, we focus on 
the asymptotic limits of the heat kernel Z(o~) = J d D xv*(x) C(x, a) and of the spectral 
dimension. For < /?* < 2, £ 2 — > when a — > 0, and the initial condition is recovered, 
since from equation (194"1) 

C(x',a) CT ~° UtorPfiv.tf)]' 1 . (122) 



Thus, Z ~ [£(a)]~ D . Still for the whole range < /3* < 2, when a — > +oo one has 
£ 2 — y +oo, $ —7- 1, and again Z ~ [£(cr)] _z? . From equations fillip and (11161) . one 
eventually gets equations (II 13[) and ( 1 1 1 7H . 

This result amends the incomplete discussion in section 14.21 With respect to 
dimensional flow, this model falls in the same class of its Hermitian-conjugate dual, 
i.e., the model with weighted Laplacian. 

j| Here we can talk about fractals because ds < da. 
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The case of g-theory with scale-dependent dimension is straightforward. One should 
simply replace £ 2 (u) in equation f 1 1 1 ft and the distributions q(x^) with a generic 
functional of, respectively, a and x M . The multi- fractional measure weight (136]) 
determines the multi-scale form of q: 



N 



q{x H 



0*00 = X^> nSgn 0Ok 



/^UM| a M 



n=l 



where we absorbed T factors into the g^nS. This also suggests the form of £ . 



N 



£ 2 (a) = Q t (a) = J^flwr* 



n=l 



The simplest case is the binomial measure f llOTft . 



(a) = k£* 



a 



a 



(123) 



(124) 



(125) 



Since < /3* < 1 in this model, at small scales aft* <C 1 the second term dominates 
and d 2 oc a 13 * , while at large scales cr/£* ^> 1 the variance of the PDF is the usual one, 



KG. 



In general, the coefficients g n may also depend on a, according to the type of profile 
one wishes to describe. Barring this possibility, from equations ( 1101ft . (1124ft . and ( 1841) . 
the spectral dimension reads 



d s (a)=D 



En WW 
J2n9n^" 



(126) 



In particular, for the case of fixed dimensionality (N = 1, 0i = (3), the spectral 
dimension coincides with equation ( 1 1 5 ft . while for the binomial measure one has 

l + A^tAx) 1 -^ 



d s (a) = D- 



(127) 



l+(4/<r)i-/3» • 

As anticipated, d$ ~ D(3* at small scales a /£* <^C 1 and d$ ~ D at large scales a/£* ^> 1. 
Since equation ( I123p is not invertible, knowing the statistics of Q does not 
immediately lead to a statistics for X. One can, however, apply equation fl 1 2 ft at 
various asymptotic regimes with different a, and give a different, patch-wise stochastic 
description of spacetime for each scale range where d$ ~ const. 



6. Discussion 

The three classes of no-scale spacetime models analyzed here display constant anomalous 
dimensions and are associated with highly non-trivial stochastic processes, all related 
(via various rescalings and mappings) to Brownian motion or scaled Brownian motion. 
Table [T] summarizes these findings. To the best of our knowledge, these processes do 
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Model 
(symmetries of C) 


d'Alembertian 


Stochastic process 


d H 


rfs 


Weighted Poincare 


eq. (J38D 


^-FSBM— V\& ) 

eq. (JZU) 


Da 
eq. ([35]) 


D(2-/3) 
eq. (PD 


Ordinary Poincare 


D., Dt = iD,(t;.) 
eqs. (|SD and (04]) 


^fsbm(c) 
eq. (EHD 


Da 


D(2-/3) 


g-Poincare 


eq. (T46D 


-^FSBM— ql^J 

eq. ||102|> 


Da 


D/3 
eq. (I105|) 



Table 1. Models of fractional spacetimes with fixed dimensionality and v = 1, 
characterized by the symmetries of the Lagrangian density, a stochastic process X(a), 
Hausdorff dimension dfj and spectral dimension ds ■ 



not have counterparts in the literature of probability theory and we had to name them 
according to the spacetime they describe. 

We also got the analytic expression of the spectral dimension of the multi-scale 
extensions of these theories, found that the models with weighted and ordinary Laplacian 
(which are Hermitian conjugate) produce the same dimensional flow, and given an 
alternative interpretation of the theory with weighted Laplacian as a fuzzy spacetime. 
In this case, dimensional flow is modified accordingly and ds — > in the UV. 

Only the model with g-Laplacian is fractal in the usual sense. Namely, the first 
Seeley-DeWitt coefficient in the heat kernel expansion corresponds to the Hausdorff 
volume and the relation between Hausdorff, spectral, and walk dimension is fl27|) . The 
other models violate these conditions, mainly because the effective density of states 
scales as the integer embedding volume. Such properties are likely to signal a difference 
in the renormalization of field theories living on these spacetimes. 

In all the models of fractional and multi-scale spacetimes we have considered in 
this paper, the spectral dimension was derived starting from statistical mechanics, via 
a generalized Langevin equation. There is, in general, an intrinsic ambiguity in the 
diffusion equation approach, inasmuch as it does not fix the the scaling dimension of 
effective diffusion time. This, in fact, is part of the definition of the theory, and one 
should not expect to get a certain value of ds without imposing a certain number of 
defining conditions. Part of this ambiguity might be related to other aspects of the 
theory which we have not considered here, such as the momentum-space structure. 
Taking as a guiding principle the conjecture jl5] according to which the spectral 
dimension ds is the dimension of momentum space, and fixing the latter by assuming 
that momentum transform is an automorphism (hence, w(k) = v(k)), one would fix the 
value of (3 a posteriori. Similar considerations are not rigorous at the present stage and 
further study will be needed. 

The above approach is classical. Quantum mechanics can provide valuable 
information about the diffusion equation, which, as a matter of fact, can be derived 
as the classical limit of a quantum diffusion process. In the appendices we discuss some 
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aspects of the relation between classical and quantum diffusion. 
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A. Quantum probability density as a bilinear 

In finding the diffusion equation (ITT]) we did not give a robust motivation to the form of 
the Laplacian K. v , i.e., of the functional form ( 170|) . We can do so by a simple calculation 
of the quantum probability density function in multi-scale spacetimes with weighted 
Poincare symmetries, closely following the same procedure as in ordinary spacetime 
[55] . Since we are in a quantum-mechanical context, we denote time with t. We also 
write D- dimensional spatial vectors in normal font. 

Consider the quantum-mechanical free particle with unit mass and Hamiltonian 
-^frcc = ~ ^sf> where /C^ is the Euclidean version of (1351) in D spatial dimensions. The 
energy E k = k 2 are given by the eigenvalue equation H iTCC e(k, x) = E k e(k, x), where we 
used the 'plane waves' ( 1391) . The associated Green equation with delta source both in 
space and time (h = 1) is 

(iV t - # frcc )G frcc (x, t; x', t') = 5 V0 (t, t') 5 v (x, x') . (128) 

The solutions are the advanced and retarded propagators 

Gl£(x,t;x>,t') = T i e fl \- *3 ] <*|e-^-'V> , (129) 

where 9 is the Heaviside step function. Transforming into energy-momentum space, and 
using 



T i6(±z)= dw^— , (130) 

it is easy to show that, 

G*£{x, t; x', t') = JdE w(E) e*{E, t)e(E, t)G$£(x, x'; E) , (131) 

where we assumed a generic measure weight w(E), we extended the 'plane waves' (1391) 
to the energy-time pair, e(E,t) = e lEt / \/vo(t)w(E) , and 

G£?(»,^fl) = /d°t m (t) ^ ( *f ■ (132) 

Making measure factors explicit and cancelling them, we obtain that the fractional 
propagators G ir ' ee are related to the ones in ordinary space G frec by 

GL A (M;*',f)= f'^m'^Zv (133) 

y/v(x)Vo{t)v(x')Vo{t') 

Comparing with the propagator computed in [T3], there is an extra prefactor l/vo(t') 
(there, the overall time-dependent prefactor is ^/v (t')/v (t)), because the propagator 
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of [H] was solution of the Green equation only with spatial delta source. Related to 
this, notice that, in the limit t — > t', G h ' cc (x, t; x', t) = 5 v (x, x')/v (t). This suggests, as a 
feedback, to define the Green equation (11281) with an ordinary time direction (ordinary 
delta and time derivative); its solutions are then \/vo(t')v (t)G h ' cc . 

Let us now turn to a generic quantum system without specifying the Hamiltonian 
H and the dispersion relation between the energy Ek and momentum. We only assume 
that H is characterized by energy eigenvalues E k and eigenstates \ipk)'- H\ipk) — Ek\i/Jk)- 
Here the momentum k can be either discrete [55] or continuous, but for simplicity we 
stick with the continuum case. Let \* x ') be a normalizable wave packet of average 
energy £ centered at the spatial point x' and with dispersion s: 



iDl „w , , n , <^-^ 



**) :=A / d D kw(k){ifj k \x')e =-|Vfc), (134) 



where A = A(x', s) is a normalization constant and we assumed a Gaussian form for 
the packet. Projecting on the position-space basis and defining ipk{%) '■= {xfyk), one 
obtains the wave-function 



(gfc-g) 2 



* x /(x) := (x\V x> ) = A / d D kw{k) ^* k (x')^ k (x)e — ^ . (135) 



Scalar products of wave-functions include a non-trivial weight factor, so that the 
orthonormality of the ipk{%) is 



(4>k>\4>k) '■= / d xv(x)ip k (x)ip* kl (x) = 5 w {k,k') . (136) 

This relation is the starting point for finding the normalization constant A by imposing 
{^x'Y&x 1 ) — 1- Defining the density of states ( 1541 . one has 

l = A 2 fdEw(E) Px ,(E)e~ iE ^, (137) 



where we made explicit the dependence of p on x' via a subscript. The result for 
A depends on the way the density of states is approximated. For instance, since 
p x '(E) oc |^(V)| 2 oc [f(x')] _1 , it is natural to replace p x >(E) with p /v(x'), where 
Po is the constant average over disorder per integer volume [55]. Thus, we obtain that 

A 2 = ^-C(s,£), (138) 

Po 

where the function C depends on the form of w(E) and, in fractional spaces, also on 

the average energy £. The details of C are not important here, since we wish to extract 

only the measure-weight dependence of the quantum PDF. 

The latter is defined as the square of the matrix element of the evolution operator 

g R (x,t;x>,t>) := ~i f^'/l A xl*-*™]**) (139) 

■\/Vo(t)Vo(t') 

d D yv(y)G R (x,t;y,t')y x ,(y), (140) 

averaged over disorder. The Hamiltonian, in fact, may be quite complicated by the 
interaction of the particle with the underlying medium. This interaction is of stochastic 
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type and may be included as an effective noise potential. We denote this stochastic 
average as ((•)). Thus, the probability density function of quantum diffusion is 

P q (x,t;x',t') := ((\g R (x,t- lX ',t')\ 2 )) 

A 2 J dEw{E)<B*(E,t)<B(E,t') ( dE'w(E')e*(E',t')e(E',t) 

(E-S) 2 + (e'-S) 2 



A 



x e 4s 

2 



dE / dE'e 



{{G R {x,x';E)G A {x',x;E'))) 

-i(E-E')(t-t')- ^~ £ ^ E '-^ 



u = E — E' 



v (t)v (f) 

x ((G R (x,x';E)G A (x',x;E'))) 

A 2 f . ,_, f , ■ ,, ,n (E-£) 2 + (,E-u,-e) 2 

dE / dwe _lw(M) e = 



v (t)v (t>) 

x ((G R {x, x'\ E)G A {x', x;E- u))) , (141) 

where we used the analogue of equation (11311) for the interacting case. At this point, we 
make two assumptions [55]: first, that u ^ s (spread much larger than the frequency) 
and, second, that 

P q {x,x';ou) := — ((G K {x,x';E)G A {x',x;E-uj))) (142) 

Po 

depends on the energy E only weakly. We can thus take the stochastic average outside 

of the integral in E. Then, equations H138[) and (11411) yield 

P q (x,t;x',t') « y } fdue-^^P^x'-u) . (143) 

Vo(t')V (t) J 

The extra prefactor l/[vo(t')v (t)} would disappear if one assumed a Green equation 
with ordinary time delta and time derivative instead of f !128[) . 

From equations (11331) . (11421) and (11431) . it follows that the coordinate dependence of 
the quantum PDF in fractional spaces is the usual one times a prefactor (spatial vector 

notation reinstated) w(x')[w(x')^(x)fo(^ / )' y o(^)]~ 1 = [ u ( x )' y o(^)]~ 1 [wo(^)]~ 1 ^ [ t; ( x )] _1 - 
The classical PDF stems from various approximations of P q [SB] , which however do not 
alter this prefactor. This is the origin of the functional form of equation (j70|) . 

B. Quantum-mechanical time 

In the previous section, we derived the PDF of a quantum diffusive process and found 
that, as in the standard case [55], it is of the form P q ~ GG, where G is the propagator 
of the particle. The classical PDF P can then be derived from P q . (It is too heuristic to 
regard the diffusion equation as the Wick-rotated version of the Schrodinger equation 
and the classical PDF P as the Euclideanized version of the quantum propagator 
G.) We thus obtained the coordinate-dependent normalization (I7D|) . This is not the 
only information one can extract from quantum mechanics. In fact, the time-space 
dependence of the propagator G, calculated from the path integral, fixes the natural 
quantum- mechanical time T, i.e., the time with respect to which the path integral yields 



CONTENTS 36 

the transition probability. Since the scaling ratio of this time with respect to spatial 
coordinates filters down to the classical level, this provides an independent identification 
of diffusion time and a check of the diffusion equations ( I7T|) and f[8"B"j) in the fixed- 
dimensionality case, where the relation between dispersion and diffusion time is given 
by equation (IB"8"|) . To this purpose, the calculation of the free-particle propagator Gf ree 
suffices. In what follows, the symbol q always denotes canonical coordinates (one spatial 
dimension for simplicity), while for geometric coordinates we reserve the symbol g(x). 

B.l. Weighted Lagrangian 

The quantum mechanics corresponding to the fractional scenarios with weighted 
Laplacians features derivatives T> in the Lagrangian. The free-particle action is 

S= [ dt"v (t")L, (144) 

where L = |m(I> t g) 2 . The Green function as a path integral 

g ~ Yl elS ( 145 ) 

paths 

was computed in [T3j for this system and it reads 



GUq ^ q ^l ^m^¥) \JM^V) exp { w^) l (146) 

The relative scaling between space and time coordinates thus roughly identifies the 
quantum-mechanical time 

which, for a power-law measure, yields 

T oc t 2 - p , (148) 

in agreement with ( |69|) in the case v = 1. Equation ( I147P is not a rigorous definition 
stemming from equation ( 114=611 . where the time-space scaling ratio is quite entangled. 
This is the reason why we do not get a precise matching with equation (1681) . which 
instead we obtain in the next example. 

B.2. Ordinary Lagrangian 

The standard-Laplacian scenario corresponds to a Lagrangian L = ^mq 2 . For this case 
we report the full calculation, which is similar to the one in [14]. Partitioning the time 
interval t — if into iV infinitesimal parts, t = t^ > tu-i > ■ ■ ■ > t\ > to = if, the action 

is 

*= IE r^vomim? ^E'M',,)%^, (149) 

n=0 ^ tn n=0 ™ +1 n 
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where in the second step one could replace v (t n+1 ) with v (t n ) or the average 
[vo(t n +i) ~ v o{tn)}/2 without changing the final result. 

To sum over all trajectories and get fll45|) . it is sufficient to integrate over all possible 
q n = q(t n ) (ordinary integration, without weights): 

N-l 



G N := K N / dgi . . . d<?jv-i exp 



y~ Y^ V o(tn+l)(qn+l ~ q n ) 



n=0 



(150) 



where K^ is a constant and we chose a partition with identical segments e = t n+ i — t r 
Ne = t — t'. It is easy to show by iteration that 

N 



Gn^n, tN'-, <?0) *o) — Kn Y 



2irie m 1 



\ y mv (t n ) V 2-iriT N 



exp 



m (qN - go) 



2i 



-JV 



where 



N 



^-E 



'm) 



(152) 



The constant Kn can be chosen so that GN(qN,qo) = / dqAiGN(qN,qM)GM(qM,qo), 
which implies Km = KmKn-m- Noticing that this relation is satisfied by 

j mv {t n ) 



K N = \{ 



?1=1 



2vrie 



we thus have 



TTl 

G N {q N ,t N ;q ,to) = \/^^ ex P 



m {q N ~ qof 



(153) 



(154) 



2i T N 
In the double limit N — > oo, e — > 0, with Ne = t — t' fixed, we finally get (q^ = q, 

qo = q') 



G bee (q,t;q',t') := lim G N (q,t;q',t' 

N— >oo 



where 



T = lim T/v 



m 1 
2^iT 

1 dt" 
, v (t")'' 



exp 



A2 



m (q — q') 
2i ¥ 



(155) 



(156) 



which coincides with equation (168J1 with v — 1. For a power-law weight Uo(t) oc t" -1 , 
equations (11561) and (J147P coincide with (11481) up to a positive constant, but otherwise 
they are different. 

B.3. q-Lagrangian 

This case is obvious, and from the ordinary Green function it follows that T = g(t). In 
the case of fixed dimensionality, T oc t 13 . 
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